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Abstract 

We propose a Bethe-Ansatz-type solution of the open spin- 1/2 integrable XXZ 
quantum spin chain with general integrable boundary terms and bulk anisotropy values 
in/{p+ 1), where p is a positive integer. All six boundary parameters are arbitrary, and 
need not satisfy any constraint. The solution is in terms of generalized T—Q equations, 
having more than one Q function. We find numerical evidence that this solution gives 
the complete set of 2^ transfer matrix eigenvalues, where N is the number of spins. 



1 Introduction 



Although the closed (periodic boundary conditions) spin-| XXZ quantum spin chain was 
solved many years ago [U 12 E] , the corresponding open chain with general integrable bound- 
ary terms has remained unsolved. Nevertheless, much progress has been achieved on this 
fundamental problem. The special case of diagonal boundary terms was solved by Gaudin |I] 
and by Alcaraz et al. Sklyanin |6 i constructed the commuting transfer matrix in terms 
of solutions of the bulk and boundary [7j Yang-Baxter equations. The general solution 
by de Vega and Gonzalez-Ruiz jH] and by Ghoshal and Zamolodchikov [9 of the boundary 
Yang-Baxter equation led directly, through Sklyanin's construction, to the general integrable 
Hamiltonian 1 



U = ft +^ sinfi?? 



coth a_ tanh j3_a\ + cosech a_ sech /3_ ( cosh 0_a x + i sinh 9_a\) 



coth a + tanh (3 + a z N + cosech a + sech (3 + ( cosh 6 + a x N + i sinh 9 + a 



XI) 



where the "bulk" Hamiltonian is given by 

JV-l 

2 



l N-l 



n=l 



Here a x , a y , a z are the standard Pauli matrices, r\ is the bulk anisotropy parameter, a>± , fi± , Q± 
are arbitrary boundary parameters, and N is the number of spins. Note that the boundary 
interactions include nondiagonal terms (proportional to a x and a y ), which originate from 
corresponding nondiagonal terms in the solutions [HI E] of the boundary Yang-Baxter equa- 
tion. Due to the presence of these nondiagonal boundary terms, the Hamiltonian does not 
have a simple reference (pseudovacuum) state. Hence, many of the standard techniques for 
solving integrable models cannot be applied to this model. 

A first step toward solving the case of nondiagonal boundary terms was taken in |lUj . 
There it was found that, for bulk anisotropy values 

171 

V = — T , p=l,2,..., (1.3) 
p+ 1 

(and hence q = e v is a root of unity, satisfying q p+1 = —1) and arbitrary values of the 
boundary parameters, the transfer matrix t(u) (see Section^ obeys a functional relation of 
order p + 1. For example, the first three functional relations are given by 

p=l: t{u) t(u + rj) - 5(u) - 5(u + rj) = f(u) , (1.4) 



1 As discussed further in Sectional here we use a parametrization of the boundary interactions which 
differs from that used in [HI E]- 
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p = 2 : t{u) t{u + 17) t(u + 2rf) - 8(u) t{u + 2rj) - <J(tt + 77) t(u) 

-5(u + 2r]) t(u + r}) = f(u), (1.5) 

p = 3 : t(u) t(u + 77) t(u + 2rj) t(u + 3rj) - 6(u + Srj) t(u + 77) t(u + 2rj) 

-5(u) t(u + 2r]) t(u + 3r?) - 5(u + 77) t(u) t(u + 3r/) (1.6) 

-5(u + 2rj) t(u) t(u + r]) + 5{u + 77) 8(u + 3rj) + 5(u) 5(u + 2r/) = f(u) , 

where 8{u) and /(w) are scalar functions which depend on the boundary parameters, whose 
explicit expressions are given in Appendix EI (Similar results had been known for closed 
RSOS models [HJII2] and closed spin chains OdlUS].) Expressions for the eigenvalues of 
the transfer matrix were also proposed in ^Oj- However, since these expressions are rather 
complicated and are in terms of zeros of the eigenvalues, rather than zeros of a Q function 
(as in conventional Bethe Ansatz), they are probably not very useful. 

By exploiting these functional relations [TH] as well as by other means [d CHI UE], a 
more conventional Bethe Ansatz solution was then found for arbitrary values of the bulk 
anisotropy, provided that the boundary parameters obey the constraint 

1 - e 2 

a_ + ei/3_ + e 2 a + + e 3 /3 + = e (6L - 9 + ) + r]k H — m mod (2m) , eie 2 e 3 = +1 , (1.7) 

where = ±1, and k is an integer such that \k\ < N — 1 and iV — 1 + k is even. A further 
drawback of this solution is that completeness is not straightforward, as two sets of Bethe 
Ansatz equations are generally needed in order to obtain all 2 N levels [18J. Nevertheless, 
finite-size effects for this model and for the boundary sine-Gordon model jU] have been 
computed on the basis of this solution [2111 • (Related results for the boundary sine-Gordon 
model have been obtained by different methods in |2I]-) Many other interesting features, 
applications and generalizations of this solution have also been found (see, e.g., 22l-[2lsj"). 

There remains the vexing problem of solving the model when the constraint (jl.7j) is not 
satisfied, i.e., for arbitrary values of the boundary parameters. Our goal has been to solve 
this problem for the root of unity case (jl.3j) . Some progress was already achieved in |2H1 ED] , 
where two of us (R.M. and R.N.) proposed Bethe- Ansatz-type solutions for special cases with 
up to two free boundary parameters (and with the remaining boundary parameters fixed to 
specific values). For those special cases (as well as for the cases where the constraint (jl.7j) 
is satisfied), the quantity A(u) defined by 

v 

A(u) = f(u) 2 -4l[6(u + jr ] ), (1.8) 

3=0 

(where f(u) and 5(u) are the functions appearing in the functional relations, e.g. (jl.4J) - 
Q1.6p ) is a perfect square. However, for generic values of boundary parameters, the quantity 
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A(u) is not a perfect square, and it had not been clear to us how to proceed. It is on this 
generic case that we focus in this paper. 

We find, for generic values of the boundary parameters, expressions for the eigenvalues 
A (it) of the transfer matrix t(u) in terms of sets of il Q functions" {cii(u) , bi(u)}, whose zeros 
are given by Bethe-Ansatz-like equations. (See (jS301), (|3~2T|) for p > 1; and (|3~T2]) . (j4~T3j) 
for p = 1.) Such generalized T — Q relations, involving more than one Q function, appeared 
already for certain special cases [SHj, and were used in [HI] to compute the corresponding 
boundary energies in the thermodynamic limit. We have verified the T — Q relations numer- 
ically for small values of p and N, and confirmed that they describe the complete set of 2^ 
eigenvalues. 

The outline of this paper is as follows. In Sectional we briefly review the transfer matrix 
and its relation to the Hamiltonian (jl.lj) . We present our solution for the case p > 1 in 
Section El and for the XX chain (p — 1) in Section HJ We conclude in Section with a 
discussion of our results. Some details are relegated to the appendix. 

2 The transfer matrix 

The transfer matrix t(u) of the model is given by |E] 

t(u) = tr Kq (u) T (u) K (u) f (u) , (2.1) 
where T (u) and T (u) are the monodromy matrices 

T (u) = R 0N (u) ■ ■ ■ R 01 (u) , f (u) = R 01 (u)---R m (u), (2.2) 
and tr denotes trace over the "auxiliary space" 0. The R matrix is given by 



R(u) 



( sinh(w + 77) \ 

sinhii sinliT? 

sinh?7 sinhw 

\ sinh(w + 77) / 



(2.3) 



where rj is the bulk anisotropy parameter; and K^{u) are 2x2 matrices whose components 
are given by [HI El 

K^y (u) = 2 (sinh a_ cosh (3_ cosh u + cosh a_ sinh (5_ sinh u) 
(u) = 2 (sinh a_ cosh (3_ cosh u — cosh a_ sinh /3_ sinh u) 
K u (u) = e e -sinh2u, (u) = e~ 9 ~ sinh 2u , (2.4) 
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and 

K+{u) = K-(-u- V )\ iaM _^ ( _ a+M) , (2.5) 

where a T , (3 T , 9 T are the boundary parameters. 2 
For u = 0, the transfer matrix is given by 

t (0) = col , Co = — 8 sinh 2Ar 77 cosh 77 sinh a_ sinh a+ cosh (3_ cosh /3 + . (2.6) 

For r] 7^ z7r/2, the Hamiltonian (jl.ljl is related to the first derivative of the transfer matrix 
at u — 0, 

ft = Cl t'(0) + c 2 I, ( 2 - 7 ) 

where 

ci = — (16 sinh 27V_1 77 cosh 77 sinh a_ sinh a + cosh j3- cosh , 

sinh 2 1] + N cosh 2 77 
2 cosh 77 

For the special case rj = in/2 (i.e., p = 1), 

t(0) = , 1/(0) = d I , d Q = (-l) N+1 8i sinh a_ sinh a + cosh /3_ cosh j3+ , (2.9) 
and the Hamiltonian (jl.lj) is related to the second derivative of the transfer matrix at u = 

ft = dit"(0), di = (-l)^ 1 (32 sinh a_ sinh a+ cosh /?_ cosh /?+) _1 . (2.10) 

In addition to the fundamental commutativity property 

[t(«),t(v)] = 0, (2.11) 
the transfer matrix also has m periodicity 

t(u + Z7r) = t(u) , (2.12) 

crossing symmetry 

t(- u - rj) = t(u) , (2.13) 

and the asymptotic behavior 

e «(27V+4)+»7(iV+2) 

t(u) ~ - cosh(#_ - 0+) for u^oo. (2.14) 



2 We use a parametrization of the boundary parameters which differs from that in |SJ|5J^]. Specifically, 
the matrices K^(u) are equal to those appearing in the second reference in |lf>| divided by the factors k^, 
respectively. 



3 The case p > 1 



We treat in this section the case (jl.3j) with p > 1, i.e., bulk anisotropy values 77 = y , y , — 
Following Bazhanov and Reshetikhin we first recast the functional relations for the 
transfer matrix eigenvalues A(u) as the condition that a matrix .M(w) have zero determinant. 
The equations for the corresponding null eigenvector, together with a key Ansatz Ij3.37|) - 
(|3.38j) . then lead to the desired set of generalized T — Q relations for A(u) (|3.4(Jj) . (|3.41j) and 
the associated Bethe- Ansatz equations ()3.43|) - ()3.50|) . 



3.1 The matrix M.{u) 

Our objective is to determine the eigenvalues A(u) of the transfer matrix t(u). As noted in 
the Introduction, the transfer matrix satisfies a functional relation (e.g., (|1.4|) - ()1.6p ). where 
the functions 5(u) and f(u) are given in Appendix [XJ By virtue of the commutativity 
property (j2.11j) . the eigenvalues satisfy the same functional relation as the corresponding 
transfer matrix, as well as the properties (j2.12j) - (j2.14|) . Hence, for example, for p = 2 the 
eigenvalues satisfy 



A(u) A (it + 77) A (it + 2rj) - 8{u) A(u + 2rj) - 5(u + rj) A(u) 

- 5(u + 2rj) A(u +rj) = f(u) . 



(3.1) 



The first main step is to reformulate the functional relation as the condition that the 
determinant of some matrix vanish. To this end, let us consider the (p + 1) x (p + 1) matrix 
M.(u) given by 



M{u) 



A(«) 
-n x (u) 



V 



-m. 



p+i 



(«) 



— mi{u) 
A(u + rj) 








-m 2 {u) 
















-n p+1 (u) 




-np-i('u) A(u+(p-l)ri) -m p (u) 
—n p (u) A{u + prj) J 



(3.2) 



where the matrix elements {rrij(u) ,nj(u)} are still to be determined. Evidently, this matrix 
is essentially tridiagonal, with nonzero elements also in the lower left and upper right corners. 
One can verify that in order to recast the functional relations as 



det.M(w) = 



(3.3) 



it is sufficient that the off-diagonal matrix elements {rrij{u) ,rij(u)} be periodic functions of 
u with period in, and satisfy the conditions 



mj(u) nj(u) = 6(u+ (j - l)rj) , j = l,2. 



,P+1 



(3.4) 
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p+1 p+1 

We now proceed to determine a set of off-diagonal matrix elements {mj(u) , rij(u)} which 
satisfies these conditions. Using (|3.4|) to express rij(u) in terms of rrij(u), and then substi- 
tuting into ()3.5|) . we immediately see that the quantity z(u) = Y\ 1 j=i m j( u ) must satisfy 

1 P 

*(«) + — - J] S(u + jrj) = f(u) . (3.6) 

*(«) , A i 

This being a quadratic equation for we readily obtain the two solutions 

zHu) = \ (f(u) ± VAR) , (3.7) 

where the discriminant A(-u) is the quantity ()1.8|) mentioned in the Introduction, 

v 

A(u) = f(u) 2 -4l[5(u + j V ). (3.8) 
j=o 

In short, we must find a set of matrix elements {rrij(u) , n,j(u)} which satisfies (|3.4j) and also 

\{m ] {u) = z ± {u), (3.9) 

where z ± (m) is given by (|3.7I) . 

In previous work P^I^9*ll3*U] we considered special cases for which A (it) is a perfect square. 
However, for generic values of the boundary parameters, A(u) is not a perfect square. Hence, 
the off-diagonal matrix elements cannot all be meromorphic functions of u. 

In order to determine these matrix elements, it is convenient to recast the expression for 
z ± (u) into a more manageable form. Noting that (see (jA.2j) . (jA.6|) and (|A.8jl ) 



we see that 



where we have defined 



i[S (u + jri) = f (u)\ (3.10) 

j=0 



A(«) = /o(«r A!(u), (3.11) 



A 1 {u) = f 1 {uf-A\{5 l {u + 3 7 1 ). (3.12) 

j=0 
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It follows from (jHHJ) and lETTTJl that 

z ± (u) = f (u)zf(u), (3.13) 

where 

= \ (/i(«) ± a/^R) • (3-14) 

Using the explicit expressions for (|A.3|I and fi(u) (|A.7jl . (|A.9J) . one can show that 
Ai(m) ()3.12j) can be expressed as 

2 

Ai(u) = 4sinh 2 (2(p+l)M)^/i fc cosh fc (2(p+l)M), (3.15) 

/c=0 

where the coefficients which depend on the boundary parameters, are given in the Ap- 
pendix flUIU), |A~T2jl for even and odd values of p, respectively. It follows from (|3.14J) and 
fl3~T5jl that 

zf{u) = X - (A(V) ± gi{u) Y(u)) , (3.16) 



where we have defined 



and 



gt{u) = 2sinh(2(p + l)u) (3.17) 



Y{u) 



\ 



^^cosh fc (2(p+ l)u), (3.18) 



k=0 



which we take to be a single-valued continuous branch obtained by introducing suitable 
branch cuts in the complex u plane. 3 One can see that Y(u) has the properties 

Y(u + r]) = Y(u) , Y(-u) = Y(u). (3.19) 

It follows from (lA~To|) . (ETTKJ) and (l3~T7j) that 

zf(u + r]) = zf{u) z+{-u) = z^{u) . (3.20) 

In short, z ± (u) is given by (|3.13|) . where zf(u) is given by ()3.16|) - (|3.18|) . and has the 
important properties (|3.20jl . 

3 We assume that the boundary parameters have generic values, and therefore, the function 
Y^ik=o ^ k cosh fe (2(p + l)u) is not a perfect square. The branch points are zeros of this function. 
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In order to construct the desired set of matrix elements, it is also convenient to introduce 
the function h(u), 

h(u) = h (u) hi (it) , (3.21) 

where h (u) is given by 



h (u) = (-If smh^(u + V )— ^— f , (3.22) 

smh{2u + r]) 

and satisfies 

ho(u) ho(-u) = S (u-rj), (3.23) 

p p 

Y[h {u + hj) = ]lh {-u-kT]) = f (u). (3.24) 

k=0 k=0 

Moreover, hi(u) is given by 4 

h^u) = (-1) N+1 4 sinh(w + a J) cosh(w + /L) sinh(w + a+) cosh(w + , (3.25) 
and satisfies 

hiu) h x {-u) = 6 1 (u-rj). (3.26) 

We are finally ready to explicitly construct the requisite matrix elements: 

mj{u) = h(-u-jri), rij{u) = h(u + jij) , j = l,2,...,p, 
z~(?z) zf (tt) h (-u) 

mp+l[u) m=i K-u - kn) ~ U.U hi(-u - k v ) ' 

z + {u) z i( u ) ho(u) 

nLi Ku + k V ) = uuhiu+kv) 

Indeed, using (|3.23|) . (|3.26|) and the fact 

p 

z + (u) z~ (u) = Yl S(u + jr]) (3.28) 

3=0 

(which follows from ()3.7|) and ([3. 8)) ). it is easy to see that the condition ()3.4j) is satisfied. It 
is also easy to see that the condition ()3.9|) (with z~(u) on the RHS) is also satisfied. We 
note here for future reference that 

n p+1 (u) = m p+1 (-u) , (3.29) 

which follows from (j3.2(J|) . We also note that if the constraint (jl.7|) with t\ = = £3 = 1 is 
satisfied, then Ylk=o^( u + krj) = z ii u )i as follows from the identity (A. 8) in Hence, 
for this case, n p+ i(u) = h(u), and the matrix M.{u) reduces to the one considered in [TT)] . 

4 Presumably, one can use the more general expression h\{u) — (— 1) +1 4 sinh(it + a_)cosh(it + 
6i/3-) sinh(u + 620!+) cosh(w + £3/?+), where = ±1, which also satisfies (13.2611 . However, for simplicity, 
we restrict to the special case e< = 1. 
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3.2 Bethe Ansatz 

The fact (13 .3j) that Jv[{u) has a zero determinant implies that it has a null eigenvector 
v (u) = (v 1 (u),v 2 (u) ,v p+1 (u)), 

M(u)v(u)=0. (3.30) 

We shall assume the periodicity 

Vj(u + in) = vj(u + (p + l)rj) = Vj(u) , j = l,...,p+l, (3.31) 

which is consistent with the periodicity M.{u + in) = M.{u). It follows from ()3.3()|) and the 
expression ()3.2|) for M.{u) that 

A(u + (j - 1)77) Vj(u) = nj-i(w) fj-i(M) + mj(w) , j = 1 , 2 , . . . ,p + 1 , (3.32) 

where Vj +p+ i = Vj and rij +p+ i = rij. Shifting u 1— ► m — (j — 1)77, we readily obtain 

A(tt) ^i(m) = h(-u - 77) t> 2 (w) + n p+1 (u) v p+l {u) , 

A(u) Vj(u - (j - l)r/) = h{u) Vj-x(u - (j - 1)t?) + - r/) v j+1 (u - (j - l)r)) , 

J = 2 , 3 , . . . , p , 

A(m) w p+ i(m — pr/) = /i(w) t> p (u — pr/) + m p+1 (M — pr/) Vi(u — pr/) . (3.33) 

The crossing properties of the eigenvalue A(— u — 77) — A(u) ()2.13|) together with ()3.29|) 
suggest a corresponding crossing property of v (u), namely, 

Vj(-u) = v p+2 -j(u) , j = 1,2,... ,p + l. (3.34) 

In particular, for j = | + 1 (which occurs only for p even !), this relation implies that v e+i(u) 
is crossing invariant, 

vv_ +1 (-u) = ve +1 {u) . (3.35) 

Moreover, (|3.34|) implies that at most |_§J + 1 components of v(u) are independent, say, 
{vi(u) , . . . , where [_ J denotes integer part. 

Substituting the explicit expression for n p+ i(u) (|3.27|) into (|3.33|) . we obtain the relations 

zf(u) h Q (u) 

nUhi(u+k V ) 

A(u) Vj(u - (j - l)//) = h(u) Vj-x{u - (J - 1)77) + h(-u - 77) v j+1 (u - (j - l)r/) 



A(u)«i(«) = - 77) 7j 2 (n) + ^ \';* K ^ v 1 {-u) , 



J = 2,...,L|j+l, (3.36) 
9 



which evidently resemble a system of generalized T — Q equations. However, since A(u) is 
an analytic function of u for finite values of u 5 , the functions Vj(u) cannot be analytic due 
to the presence of the z1[(u) factor in ()3.36|) . 

We therefore propose instead the following Ansatz: 

Vj (u) = a 3 (u) + b,{u) Y{u) , j = 1 , 2 , . . . , ||j + 1 , (3.37) 
where Y{u) is the function (|3.18j) . and cij(u) , bj(u) are periodic, analytic functions of u, 



CLj [u ) 



Aj J] smh(w - , bj(u) = B j J] sinh(u - «<*>) , j ^ V - + 1 



2M a 2M 6 

[ S * 

1 

n 

fc=] 

Mb 



k=l k=l 

M a 



sinh(u — u k ^ ) sinh(u + u k 



= sinhtw — «[, 7+1 ) sinh(w + u k 7+1 ) , (3.38) 



k=l 



whose zeros {u k a ^ , uf 1 ^}, normalization constants {Aj , Bj}, and also the integers M a , are 
still to be determined. 6 The forms (I3.38J) for a,j(u) and bj(u) evidently have the periodicity 
and crossing properties 

a,j(u + m) = dj{u) , bj{u + in) = bj{u) , j = 1 , . . . ,p + 1 , 

0|+i(-u) = a|+i(u), 6f+i(-w) = 6f+i(w) , (3-39) 

which reflect the corresponding properties of Vj(u) ()3.31|) . ()3.35|) and of Y(u) ()3.19|) . We 
have obtained numerical support for this Ansatz, which we discuss at the end of this section. 

We now substitute the Ansatz ([3.37)1 . as well as the expression for zf(u) ()3.16jh into 
(I3.36j) . Since A(u) and Y{u) 2 (but not Y(u) !) are analytic function of u, we can separately 
equate the terms that are linear in Y(u), and the terms with even (i.e., or 2) powers of 
Y(u). In this way we finally arrive at the generalized T — Q equations: 

A(u) ai(u) = h(-u - rj) a 2 (u) + p \fi(u) ai(-«) + F(u) 2 &i(—u) 

z llk=i n n u + Kr l) L 

A(u) %(it - (j - l)r/) = h{u) Oj_i(m - (j - 1)77) + /i(-tt - 77) a J+ i(w - (j - 1)77) , 

J = 2,... +1, (3.40) 
5 This is a well-known consequence of the transfer matrix properties (|2.11ll - (|2.14|l . 

6 Since the normalization of the null eigenvector v(u) is arbitrary, one of the normalization constants, say 
Bi, can be set to unity. 
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and 

A(u) bi (u) = h(-u - 7]) b 2 (u) + fttu) h(-u) + g-i(u) oi(-m) 

A(u) - (j - 1)77) = /i(u) 6j_i(u - (j - 1)77) + - 77) 6 i+ i(M - (7 - 1)77) , 

J = 2,...,L|j+l, (3.41) 

where cle +2 (u) = az(—u) and a,p+3 (w) = a P +i (—u) for even and odd values of p, respectively, 
and similarly for the b's. 

The asymptotic behavior A(u) ~ e «( 2Ar + 4 ) f or M __>. go (|2,14j) together with the T — Q 
equations imply the relation 

M a = M b + p+l. (3.42) 

An analysis of the ti-independent terms yields relations among the normalization constants 
and sums of zeros (Yli u i > Y^i u i)i which we do not record here. 

As usual, analyticity of A (it) and the T — Q equations imply Bethe-Ansatz-like equations 
for the zeros } of the functions {aj(u)}, 

ho(-u^ - rj) = h{uf l] ) aii-vth +9i(u ( ? l) ) Yjv^T hj-u^) 
h (u^) 2a 2 {u^) hi-u^ - V ) UU hM ai) + kri) ' 

h(-u[ aj) - jrj) a i _i(u, (a>) ) p 



j = 2,... ,£J+1, (3.43) 



%^ ) + (i-l)77) """' L 2- 
and for the zeros {itf* } of the functions {bj(u)}, 



h (u^) 2b 2 (u^) hi-u^ - 77) nLi hi{u? l] + kq) ' 

K-^-M - 6 i- 1 ^ ) ) J = 2,...,^j+1. (3.44) 



%; 6f) + a - 1)77) ' ' 2- 

Moreover, there are additional Bethe-Ansatz-like equations for the normalization constants. 
Indeed, noting that ho(u) has a pole at w = — ?, it follows from the analyticity of A (it) and 
the T — Q equations ()3.40|) that 

«i(|) = a 2 (-|), (3.45) 

= J = 2,-.. +1. (3.46) 

11 



In obtaining the first equation (|3.45jl . we have made use of the identity 

/i(-|)=2n/»i(-|+^). 

k=0 



(3.47) 



The equations (|3.45|1 . ()3.46|) evidently further relate the normalization constants {A/}- Sim- 
ilarly, the T — Q equations (j3.41|) imply 



61 ( 



Lfj + 1 



(3.48) 
(3.49) 



which relate the normalization constants {Bj}. Finally, noting that the first (i.e., j = 1) 
T — Q equation in the set (J3.41)) has the factor YH=i h\{u + krj) in the denominator which 
can vanish, e.g. at u = — cc_ — 77, leads to the relation 



- rj) 61 (a_ +rj) = -</i(-a_ - 77) oi(a_ + 77) 



(3.50) 



which relates the normalization constants A\ and B\. A similar analysis of the first equation 
in ()3.40|) gives an equivalent result, by virtue of the identity /i(«o) 2 = gi(uo) 2 Y (uq) 2 if uo 
satisfies 11^=0^ ( w o + JV) = 0> which follows from ()3.12|1 and the fact (j3.15|) that Ai(u) = 
g x {ufY{uf. 

The energy eigenvalues of the Hamiltonian (jl.lj) follow from (|2.fij) - (|2.8j) and the T — Q 
relations IpOITjl . 



£ = CiA'(O) + c 2 = C!C 



a;-(-(j - l) V ) + a^iHj ~ 1)^7) + /i'(0) 



+ c 2 



(3.51) 



a,(-(j - 1)77) a^i-U - l)r)) h(0) 
where j can take any value in the set {2 , . . . , |_§J + 1}. For j ^ | + 1, it follows that 

2M a 

E = - sinliT? ^ [coth(w[ %) + (j - 1)77) - coth^ -0 + (j - 1)77) 
z=i 

1 1 
+ - sinhr/ (cotha_ + tanh/3_ + cotha + + tanh/9 + ) + -(A — 1) cosh?7 , (3.52) 

which does not depend explicitly on the normalization constants. For j = | + 1, there is 
an additional contribution from the term aj [ £ ]\ v l in (|3.51|) . since this cij(u) is crossing 



invariant. It follows that 

U 

sinh n < 

2 



* 1 x. 

E = -smh?7 <^ 2^ 



=1 



«j(-C/--l)»0 



coth(4 + y) - coth(uJ - y; 



2iU„ 



1=1 



coth( M , * + y ) |> + 



(3.53) 
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where the ellipsis denotes the terms in (I3.52|) that are independent of Bethe roots. If one 
works instead with the T — Q relations ()3.41|) . one obtains the same results ()3.53|) . ()3.52j) . 
except with sums over the b roots. 

We have verified the T — Q equations numerically, for values of p and iV up to 6 and for 
generic values of the boundary parameters, along the lines [T%] . These results are consistent 
with the conjecture 

M a =[^-±\+2p + l, M 6 =L^^J+P, (3.54) 

which agrees with the relation ()3.42|) . These numerical results also indicate that our Bethe 
Ansatz solution is complete: for each value of N, we find sets of Bethe roots corresponding 
to each of the 2 N eigenvalues of the transfer matrix. As already remarked, this numerical 
work provides support for the Ansatz ()3.37|) . (j3.38|) . 



4 The XX chain (p = 1) 

The case p = 1 corresponds to bulk anisotropy value rj = in/ '2, for which the bulk Hamilto- 
nian (jl.2j) reduces to 

Z n=l 

which is known as the XX chain. The open XX chain with nondiagonal boundary terms was 
studied earlier in 



The functional equation for the case p — 1 is given by (jl.4|) . We find that a suitable 
matrix M{u) is given by 

M{n)=\ A( "\ :™ lW I , (4.2) 



— rii(u) A(u + rj) 



where 



wi(m) = -7—7- — r \h (u) 5 x {u + rj) + h (-u-r]) z 1 (u)] , (4.3) 

ni(u) = \h (-u) 5t(u + r]) + h (u + rj) zf(u)] . (4.4) 

hi(u) L J 

Indeed, one can verify that the condition det J\A (u) = reproduces the functional equation 
dUJ). Note that 

n\{u) = rni(— u) . (4.5) 
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The corresponding null eigenvector v(u) = (vi(u) ,v 2 (u)) satisfies Ji4(u) v(u) = 0, i.e., 

A(u) vi{u) = mi(u) v 2 (u) , (4.6) 

A(«) v 2 (u — rj) = rii {u — rj) v\ {u — rj) . (4.7) 

The crossing symmetry A(— u — rj) — A(u) and ()4.5|) suggest 

V2 (u) = v x {-u) . (4.8) 

That is, only one component is independent, say, v\{u). Substituting the explicit expression 
(J4.3|) into the first equation ()4.6|) . we obtain 

1 



A(-u) vi(u) 



[h (u) 5i(u + r]) + h (-u - rj) z l (u)] vi(-u) . 



hi(—u 

Similarly to the p > 1 case, we make the Ansatz 

vi(u) = ai(u) + bi(u) Y(u) , 

where 

2M a 2M b 

a\{u) = A J^J sinh(-u — ) , bi(u) = sinh(-u — uf ) . 

k=l k=l 



(4.9) 



(4.10) 



(4.11) 



Substituting this Ansatz, together with the expression for z 1 {u) (|3.16|) . into (|4.9j) . we obtain 
the desired generalized T — Q equations: 



A(u) ai{u) hi(—u) 



A(u) bi(u) hi(-u) 



h (u) 5i(u + r]) + ^h (-u - rj) fi(u) 



-h (-u - rj) g x {u) Y(u) 2 h{-u) 



h (u) + rj) + -h (-u - rj) fi(u) 



ai(-u) 



(4.12) 



h(-u) 



1 



- ^ h o(-u - rj) g x {u) ai(-u) . 
From the asymptotic behavior we obtain the relation 

M a = M b + 2 . 
The corresponding Bethe Ansatz equations are 



h (-u 



(ai) 



rj) 



+ V ) ai (-«r ; ) 



(ai)> 



h (u^) 



h(u^) ai (-ut ] ) - 9i(ut ] ) h(-u^) 
25 l (u^ + V ) hi-uf^) 



a(u^) b^-un - giiuD ai (- u r) 
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Xblh 



Xblh 



,0nh 



(4.13) 
(4.14) 

(4.15) 
(4.16) 



and, e.g., 



/l(a-) h(-a-) = gi(a-) ai(-a_) . (4.17) 

There are no additional relations arising from analyticity at u = — Q analogous to ()3.45|) . 
f^OHjl due to the identity = 2<5 1 (|). 

The energy eigenvalues of the Hamiltonian (jl.l|) follow from (J2.9j) . (|2.1U|) and the first 
T -Q relation (l4~T2l . 

a[(0) , fci(O)" 



£ = diA"(0) = did 

mth 7/i ai) + 

2 



'ai(O) /n(0). 



% S y^ j cothw[ ai ^ + - (coth«_ + tanh /3_ + cotha + + tanh . (4.18) 



Working instead with the second T — Q relation ()4.13|) gives the same result ()4.18|) except 
with sums over b roots. 

We have verified that the above T — Q equations are well-satisfied numerically, for values 
of N up to 8 and for generic values of the boundary parameters, along the lines ^H] • These 
results are consistent with the conjecture 

N - 1 N - 1 

M a =L^-J+3, M fe =L^-J+l, (4.19) 

which agrees with the relation f)4. 14j) , and in fact also with (|3.54j) . These results also indicate 
that our Bethe Ansatz solution is complete. 



5 Discussion 

We have found a Bethe- Ansatz-type solution of the open spin- 1/2 integrable XXZ quantum 
spin chain with general integrable boundary terms at roots of unity. All six boundary param- 
eters are arbitrary. In particular, the boundary parameters need not satisfy the constraint 
f)1.7|) that arose in previous work |H)]-PH]. Moreover, in contrast to that earlier solution, our 
new solution appears to give the complete set of 2 N eigenvalues in a straightforward manner. 
This solution is essentially the same for both even and odd values of p, the main difference 
being that, in the former case, one of the Q functions is crossing invariant. 

Part of the price paid for this success is that there are multiple Q functions {a,j(u) , bj(u)} 
and corresponding multiple sets of Bethe roots {wj° J \ uf 1 }. However, we have already 
demonstrated the feasibility of performing thermodynamic (N — > oo) computations with two 
such sets of Bethe roots [HI] • Hence, we expect that this multiplicity of sets of Bethe roots 
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will not cause significant computational difficulty. A further complication is the appearance 
of normalization constants {Aj ,Bj} and their corresponding Bethe-Ansatz-type equations. 

Another part of the price paid for this success is that the bulk anisotropy parameter 
is restricted to the values ()1.3|) . However, we expect that it should be possible to further 
generalize our solution to the case rj = mp' /(p + 1), where p' is also an integer. Indeed, 
we expect that functional relations of order p+1 with the same structure (e.g., (jl.4j) - (jl.6|) ) 
will continue to hold for that case, except with a different function f(u) that will now 
depend also on p' . Hence, to the extent that a number can be approximated by a rational 
number, this approach should in principle solve the problem for general imaginary values 
of rj. Unfortunately, this approach does not seem to be suitable for directly addressing the 
problem of real values of rj, for which case the transfer matrix presumably does not obey 
functional relations of finite order. Nevertheless, as in the case of the sinh-Gordon and sine- 
Gordon models, it may perhaps be possible to obtain results for real values of rj from those 
of imaginary values of rj by some sort of analytic continuation. 

Although we have considered here the case of generic values of the boundary parameters 
for which the quantity A(u) (jl.8j) is not a perfect square, we find numerical evidence that 
our solution remains valid when A(u) becomes a perfect square. Presumably, for such 
special cases, the Q functions {ctj(u) ,bj(u)} are not independent. It may be interesting to 
determine the precise relationship between these Q functions and those appearing in the 
previously found solutions [23 ED]- 

We remark that the set of off-diagonal elements ()3.27|) of the matrix J\4 (u) is not unique. 
Indeed, we have found other sets of matrix elements which also give det M.{u) = 0. Among 
all the sets which we found, the particular set presented here has several advantages: (i) it 
works for both even and odd values of p; (ii) the corresponding T — Q relations and Bethe 
Ansatz equations are relatively simple; (iii) the corresponding values of M a and M& are 
minimized. Nevertheless, it may be worthwhile to continue looking for alternative sets of 
off-diagonal matrix elements, which may further reduce the values of M a and M&, or which 
may have other nice properties. 

A key step in our analysis is the Ansatz ()3.37j) . (|3.38j) . which allows us to express the 
non-analytic quantities {vj(u)} in terms of analytic ones {a,j(u) ,bj(u)}. We have numerical 
evidence that this Ansatz is valid. However, it is not clear whether this Ansatz is the most 
"economical": there may be alternative Ansatze which introduce fewer Q functions. For 
example, there may be some fixed relation between a,j(u) and bj(u). 

The structure of our generalized T — Q equations bears some resemblance to that of the 
conventional TBA equations of the XXZ chain Presumably, this common structure has 
its origin in the fusion rules and root of unity properties of the underlying U q {su2) algebra. 
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Having in hand an exact solution of a model with so many free boundary parameters, 
one can hope to be able to analyze a plethora of interesting boundary behavior. 

Finally, we note that it should be possible to generalize the approach presented here to 
open integrable anisotropic spin chains constructed from R and K matrices (both trigono- 
metric and elliptic) corresponding to higher- dimensional representations and/or higher-rank 
algebras. 

We hope to be able to address some of these issues in future publications. 
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A Appendix 

Here we collect explicit expressions for various quantities appearing in the text. 

The functions S(u) and f(u) are given in terms of the boundary parameters a T , (3 T , 6 T 



by |22| 



S(u) = 5 (u)Si{u) , 



(A.l) 



where 



S (u) 




(A.2) 



Si(u) 



2 4 sinh(u + r\ + a J) sinh(w + r\ — a) cosh(u + rj + f3J) cosh(u + rj — (3-) 



x sinh(-u + 7] + a+) sinhf/u + 77 — a+) cosh(w + r] + (3+) cosh(-u + r] — (3 + ) , (A. 3) 



and therefore, 



8j(u + in) = Sj(u) 



5j(-u-2r)) = 6j(u), 



J = 0,1. 



(A.4) 



Moreover, 



/(«) = /o(«)/i(«) ■ 



(A.5) 
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For even values of p, 

/„(«) = (-l) N+1 2-^smh w {(p+i)u) t (A.6) 
h(u) = (-l)^ +1 2 3 -^( 

sinh ((p + l)a_) cosh ((p + sinh ((p + l)a + ) cosh ((p + cosh 2 ((p + l)u) 

- cosh ((p + l)a_) sinh ((p + cosh ((p + l)a+) sinh ((p + 1)/?+) sinh 2 ((p + l)tt) 

- (-1)* cosh ((p + 1)(0_ -0+)) sinh 2 ((p+ cosh 2 ((p+ l)u) ) . (A.7) 
For odd values of p, 

/ ( u ) = (-l) iV+1 2- 2pJv sinh 2iV ((p+l) M )tanh 2 ((p+l)M) , (A.8) 

A(tt) = -2 3 ^( 

cosh ((p + l)a_) cosh ((p + cosh ((p + l)aq_) cosh ((p + l)/3 + ) sinh 2 ((p + l)w) 

- sinh ((p + l)a_) sinh ((p + l)/3_) sinh ((p + l)aq_) sinh ((p + l)/3 + ) cosh 2 ((p + l)w) 
+ (-1)* cosh ((p+l)(0_-0 + )) sinh 2 ((p+l)u) cosh 2 ((p + l)u) ) . (A.9) 

For both even and odd values of p, these functions have the properties 

f J (u + r ] ) = f j (u), f j (- u )=f j (u), j = 0,1. (A.10) 

The coefficients //& appearing in the function ()3.18|) are given as follows. For p even, 

= 2- 4p j-l-cosh 2 ((p + l)(#_-# + )) 

- cosh(2(p + l)a_) cosh(2(p + l)a + ) + cosh(2(p + l)a_) cosh(2(p + 
+ cosh(2(p + l)a + ) cosh(2(p + l)/3_) + cosh(2(p + l)a_) cosh(2(p + l)/3 + ) 
+ cosh(2(p + l)a+) cosh(2(p + 1)(3 + ) - cosh(2(p + l)/3_) cosh(2(p + 1)/?+) 
+ cosh((p + l)(a_ + a+)) cosh((p + l)(/3_ - /?+)) 

1 2 

- cosh((p + l)(a_ - a+)) cosh((p + + (3+)) 

+ 2(-l) 7V cosh((p + l)(0_ -0+))[cosh((p+l)(a_ - a+)) cosh((p + l)(/3_ -/?+)) 

- cosh((p + l)(a_ + «+)) cosh((p + l)(/3_ + /?+)) 



/'o 



//! = 2 1 ~ 4p j cosh((p+ l)(a_ - a+)) cosh((p + l)(a_ + a+)) 



+ (-l) Jv cosh((p + + /3 + ))cosh((p+ - 0+))] 

- cosh((p + - [cosh((p + + 

+ (-1)* cosh((p + l)(a_ + «+)) cosh((p + 1)(0_ - 0+)) 

Ht = 2- 4p sinh 2 ((p+l)(^-0 + )). 



(A.ll) 



For p odd, 

/i = 2^|-l-cosh 2 ((j9+l)(^-0 + )) 

- cosh(2(p + l)a_) cosh(2(p + l)a+) - cosh(2(p + l)a_) cosh(2(p + 

- cosh(2(p + l)a+) cosh(2(p + - cosh(2(p + l)a_) cosh(2(p + 

- cosh(2(p + l)a+) cosh(2(p + l)(3 + ) - cosh(2(p + cosh(2(p + l)/3+) 
+ cosh((p + l)(a_ + a+)) cosh((p + - /?+)) 

1 2 

+ cosh((p + l)(a_ - «+)) cosh((p + + /?+)) 

- 2(-l) Ar cosh((p + l)(0_ -0+))[cosh((p+l)(a_ -a+))cosh((p + l)(/3_ - 
+ cosh((p + !)(«_ + «+)) cosh((p + + /?+)) 



/ii = 2^ 4p | cosh((p+ l)(a_ + £*+)) cosh((p+ l)(a_ - «+)) 

+ (-1)" cosh((p + - cosh((p + - 

+ cosh((p + + (3 + )) [cosh((p + - 

+ (-1)* cosh((p + l)(a_ - «+)) cosh((p + 1)(0_ - 0+)) 
/i 2 = 2- 4p sinh 2 ((j9+l)(^-0 + )). 
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